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Abstract 

Cavity optomechanical systems can be used for sensitive de¬ 
tection of mechanical motion and to control mechanical res¬ 
onators, down to the quantum level. The strength with which 
optical and mechanical degrees of freedom interact is defined by 
the photon-phonon coupling rate go, which is especially large 
in nanoscale systems. Here, we demonstrate an optomechanical 
system based on a sliced photonic crystal nanobeam, that com¬ 
bines subwavelength optical confinement with a low-mass me¬ 
chanical mode. Analyzing the transduction of motion and effects 
of radiation pressure we find a coupling rate go/27r ~ 11.5 MHz, 
exceeding previously reported values by an order of magnitude. 
Using this interaction we detect the resonator’s motion with 
a noise imprecision below that at the standard quantum limit, 
even though the system has optical and mechanical quality fac¬ 
tors smaller than 10^. The broad bandwidth is useful for appli¬ 
cation in miniature sensors, and for measurement-based control 
of the resonator’s motional state. 

The motion of a mechanical resonator can be read out with extreme sen¬ 
sitivity in a suitably engineered system whose optical response is affected by 
the displacement of the resonator. The resultant coupling between optical and 
mechanical degrees of freedom also gives rise to a radiation pressure force that 
enables actuation, tuning, damping, and amplification of the resonator, with 
applications ranging from classical information processing to quantum control 
of macroscopic objectSH Such control can be established either passively, by 
employing the intrinsic dynamics of the systenPH^, or actively, by using the out¬ 
come of displacement measurement^. Fast, sensitive measurement of nanome¬ 
chanical displacement can as such be used for optic al coo linjUll, squeezed light 
generatioi P, qu antum non-demolition measurementJHE3 and enhancing sensor 
bandwidtbEIEH 

In a cavity optomechanical system, which has an optical resonance frequency 
ujc that depends on the position of a resonator, both the sensitivity of a displace¬ 
ment measurement and the magnitude of effects caused by radiation pressure 
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forces are governed by two parameters: on the one hand the strength with which 
acoustic and optical degrees of freedom interact, expressed as the magnitude of 
the resonator’s influence on the frequency and on the other hand the cavity 
linewidth n. The interaction strength is characterized at the most fundamental 
level by the vacuum optomechanical coupling rate as it enters the optome¬ 
chanical interaction Hamiltonian ^int = Tig{)d)a(h + 6’*'), where a and h are the 
photon and phonon annihilation operators, respectively. As this Hamiltonian 
shows, ^0 is the frequency response of the optical cavity due to the mechanical 
displacement in a typical quantum state, where the total number of phonons is 
of the order of 1. 

Per photon in the cavity, the effective optomechanical measurement rat^lEl, 
as well as the radiation-pressure induced alteration of a resonator’s frequency 
and damping through dynamical backaction, scale with g^/n. Improving this 
ratio is thus desirable for more sensitive measurements and for better optical 
control of the mechanical resonator. Decreasing the optical damping to a low 
value has been very fruitful, but can present several drawbacks as well: narrow 
linewidths place stringent demands on excitation sources and fabrication tol¬ 
erances, and make integration of many devices, e.g. in practical sensor arrays, 
difficult. Moreover, dynamical instabilities and nonlinear linewidth broadening 
limit the number of photons with which a high-Q cavity can be populated. Fi¬ 
nally, several schemes for me asurement and control in fact rely on fast, broad¬ 
band optical responsJlElEH The photon-phonon coupling rate g^^ vice versa, 
is given by go = Gxzpf, where G = du^jdx is the frequency shift per unit 
displacement x and = yjTip2xri^l^ are the zero-point fluctuations of a 
resonator with mass rUeff and frequency The magnitude of go is maximized 
in suitably engineered miniature systems, as G and Xzpf benefit from small cav¬ 
ity size and small resonator mass, respectively. Indeed, the highest values of go 
to dat e have been achieved in m icrom eter-size devices such as photonic crystal 
cavitie^MlMUl or disk resonator £3221, with reported values ranging up to about 

^/o/27r 1 MHJ22I2I1 

In this work, we show that optomechanical coupling rates can be signifi¬ 
cantly enhanced by using photonic modes with subwavelength confinement. We 
realize a sliced photonic crystal nanobeam in which light is highly confined in a 
nanoscale volume near the moving dielectric interfaces of a low-mass resonator, 
leading to unprecedented interaction strengths. We use a simple free-space opti¬ 
cal setup to address the structure and demonstrate optical tuning of the mechan¬ 
ical resonance frequency, as well as sensitive readout of mechanical motion. The 
observed optical forces and measurement sensitivity provide ns with two inde¬ 
pendent ways to determine the vacuum coupling rate to be ^o/27r ~ 11.5 MHz. 
We demonstrate displacement readout with a detection imprecision below that 
at the standard quantum limit, i.e. with a noise level that is comparable to 
the quantum fluctuations of the resonator. We achieve this using only 22 pW 
of detected power even in a system with modest optical and mechanical qual¬ 
ity factors. The operation with a relatively large cavity bandwidth is especially 
attractive for system integration and miniature sensor technologies as well as 
measurement-based control in nano-optomechanical systems. 

Working principle. To realize a large photon-phonon coupling rate go = 
Gxzpf, we develop a novel system that is based on a patterned silicon photonic 
crystal nanobeam, which combines optical confinement with flexural mechanical 
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motion (Fig. [^. The beam is ‘sliced’ through the middle such that it mechani¬ 
cally resembles a pair of doubly clamped beams, coupled through the clamping 
points at the ends of the nanobeam. Figure shows the simulated fundamen¬ 
tal in-plane mechanical resonance of the sliced nanobeam structure. The small 
width (80 nm) of the narrowest parts of the half-beams ensures both the mass 
(^ 2.4 pg) and the spring constant of the nanobeam are small, leading to large 
zero-point fluctuations Xzpf. 

The in-plane mechanical motion causes a strong change in the separation 
distance d. This leads to a giant change of the optical response, as we design the 
beam to concentrate light in the subwavelength slit separating the two halves. 
In general, a displacement-induced frequency shift of an optical mode depends 
on the fraction of the energy density that is located near the moving dielectric 
boundarie^^. To maximize this effect, we rely on the high localization of energy 
that can occur in systems with dielectric discontinuities with subwavelength 
dimension^, in this case provided by the narrow slit through the middle of 
the sliced nanobeam. The periodic patterning of the beam creates a photonic 
crystal, with a quasi-bandgap for TE-polarized modes guided by the beam (see 
Supplementary Information). The waveguide mode at the lower edge of the 
band gap has strongly confined electric fields in the nanoscale gap separating 
the ‘teeth’ of the two half-beams (Fig. [^). 

The truly subwavelength character of this waveguide is revealed by calculat¬ 
ing its effective mode area, which we suitably define as A = f dV W (r)/(alFniax), 
where the energy density W{r) = e(r)|E(r)p has its maximum IFmax just at 
the vacuum side of the gap boundary, and we integrate over a full unit cell with 
period a. The mode area is only 2.38 x 10“^^ m^ for a gap width of 60 nm, 
or in other words A = 0.0llA^, with A the wavelength in vacuum. In fact, it 
is even 8 times smaller than the squared wavelength in silicon^ even though 
the maximum energy density is actually localized in the vacuum gap (Fig. &)• 
This subwavelength mode area is essential to the sliced nanobeam and makes 
it stand out with resp ect to other designs, including the related double-beam 
‘zipper’ cavit 3 EElElEll, where the optical cavity modes of two photonic-crystal 
nanobeams are coupled by placing the beams close together. 

Recently it was shown that with a similar approach photonic crystal nanobeam 
cavitie s can be created that have a high quality factor and an ultrasmall mode 
volum^mH We introduce a defect in the periodicity in the middle of the beam 
so confined cavity modes are created with a frequency in the bandgap. These 
are derived from the band of interest (i.e. the lower bandgap edge) by reducing 
the width of the central pair of teeth, such that the effective refractive index is 
locally reduced. Figure shows the simulated field profile of the lowest-order 
optical cavity mode. 

Numerical simulations confirm that the frequency of both the band edge 
and the defect cavity mode derived from it respond strongly to a displacement 
of the two half-beams, reaching G = du^jdx ^ 27r x 0.4 THz/nm for a gap 
width of 60 nm (Fig. Et)- As expected, this value increases for sm aller gap sizes, 
due to an increase of the fraction of the energy in the gapESEH We define the 
displacement coordinate as x = d/2, such that it can be directly related to 
the maximum lab-frame displacement of the antisymmetric mechanical mode 
depicted in Fig. Et- Note that the choice of the definition of x is in principle 
arbitrary (with a properly matched definition of rueff), whereas the coupling rate 
^0 is independent of this definition. To determine the optical frequency shift, the 
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Figure 1 | Geometry and resonances of the sliced nanobeam, (a) Simulated 
displacement profile of the fundamental (in-plane) mechanical resonance of the struc¬ 
ture. (b) Cross section in the center of the sliced nanobeam (indicated by the dashed 
line in c) , showing the simulated energy density distribution of the fundamental opti¬ 
cal cavity mode of the structure, (c) Simulated transverse electric field profile of the 
fundamental optical cavity mode of the structure, (d) Simulated frequency shift as a 
result of an outward displacement of 1 nm. The cavity mode shift was determined by 
simulating the full nanobeam and introducing a uniform displacement along the beam, 
(e) Electron micrograph of a fabricated device. 


entire half-beams are displaced in the simulation. The displacement of the actual 
mechanical mode is not uniform along the beam (Fig. [^), meaning that due 
to the finite extent of the optical mode the value of G will be slightly reduced. 
Taking into account the optical and mechanical mode profiles (Figs. [^, 0)5 we 
estimate it to be 0.90 times the value shown in Fig. [T (see Supplementary 
Information). 

Using standard lithography techniques (see Methods for details), we realize 
sliced nanobeams in silicon with a length of 11 pm separated by an average gap 
size of 60 nm. An electron micrograph of a fabricated device is shown in Fig.[^. 

Pree-space readout. We address our structure using a simple reflection 
measurement, schematically shown in Fig. Ilf- The employed resonant scattering 
techniqu^^ places the sample between crossed polarizers to allow the detection 
of light scattered by the cavity mode (whose dominant polarization is oriented 
at 45° to the polarizers) while suppressing light reflected by the substrate. By 
scanning the frequency of a narrowband laser we record the reflection spectrum, 
depicted for one of the samples in Fig. If- The dispersive lineshape is caused by 
interference of the resonant scattering of the cavity with non-resonant scattering 
by the n anobe am. The cross-polarized reflectance R is thus well fitted by a Fano 
lineshapJHEHl: 


R{A) = 


ce 


i(f _ 


\/^in^out 


—iA + k/2 


( 1 ) 


where c and cp are the amplitude and phase of the non-resonant scattering, 
respectively, and A = cj — cjc is the detuning of the laser frequency uj from the 
cavity resonance with linewidth k. The rate at which light can couple to the 
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cavity mode from the free-space input beam is given by /^in, whereas /^out is 
the rate at which the cavity decays to the radiation channels that are detected 
through the output analyzer. In principle, these coupling rates can be unequal 
because the light emitted by the cavity has a spatial mode profile that differs 
from the Gaussian input beam. 

Fitting equation 0 to the reflection spectrum yields the center frequency 
and linewidth of the cavity, as well as a value for We determine 

/^ex ~ 0 . 29 k and the optical quality factor Qopt = oJc/k ^ 400. The measured 
Qopt is 2-3 times lower than the simulated one, a discrepancy that we attribute 
to fabrication imperfections. 

Thermal motion of the nanobeam Sx modulates the cavity frequency by 
SuJc = {duJc/dx)Sx. This produces a change in detected power proportional 
to the derivative of the reflection spectrum: 6P = Piri{dR/duJc)^^c- Here we 
assumed the intracavity amplitude is instantaneously affected by the mechanical 
motion, which is justified since ftrn (see Supplementary Information for the 
more general case). Thus, the power spectral densities of x and P are related 
as 


(I!)' (i)' <^> 

Figure 0 shows the detected spectral density Spp for the laser tuned to 
the optical resonance frequency, with a relatively high optical power incident 
on the sample (Pin = 367 pW), corresponding to a detected power of 22 pW. 
Because of the linear relation between Spp and Sxx shown in equation ([^, this 
is a direct measurement of the spectrum of thermal motion in the nanobeam. 

The two peaks at 2.6 and 3.2 MHz correspond to the two fundamental in¬ 
plane modes of the coupled halves of the nanobeam. For a perfectly symmetric 
structure, the coupling leads to two eigenmodes: a common mode, for which the 
half-beams move in phase and their separation d is not affected, and a differential 
mode, for which anti-phase movement of the half-beams results in maximal 
variation of d. Fabrication-related imperfections can break the symmetry of 
the system, such that the actual normal modes are linear combinations 
of the half-beam eigenmodes ' 0 a = HQ,(' 0 isin^ — 02 cos and = 

74 ^( 0 icos 6 > + 02 sin 6 >), where 0 can in principle take any value. As we show 
in the Supplementary Information, the splitting between the mode frequencies 
is enhanced due to the presence of compressive stress in the studied sample, 
which also reduces the mode frequencies with respect to the simulated value in 
absence of stress of 6 MHz. Since the two modes generally affect the separation 
d differently, they have different photon-phonon coupling rates which are 
maximal for a purely differential mode ( 6 > = 7 r/ 4 ). With our definition x = d/2, 
this is reflected in the fact that the ratio between the zero-point fluctuation 
amplitudes of the normal modes is = y/f^/ 3 (l + sin 2 d)/r^Q,(l — sin 2 d) 

(see Supplementary Information). The variance in x due to thermal motion in 
the two modes is set by the equipartition theorem, taking into account this 
difference in Xzpf. The ratio between the areas of the two resonance peaks in the 
experimental spectrum of Spp therefore directly yields the mixing angle 0. 

In fact, fitting two resonant modes to the displacement spectrum also allows 
determining the transduction factor that relates the measured optical power 
spectral density Spp io the displacement spectrum Sxx- To do so, we calculate 
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the thermal variance (x^)th = We determine from the 

measured 0 and from the effective mass of purely antisymmetric motion, which 
we computed from the simulated displacement profile to be rueff ~ 0.39m, with 
m the total mass of the beam. We further assume that the temperature T of the 
mechanical bath is equal to the lab temperature. The validity of this assumption 
is tested by performing power- and detuning-dependent measurements presented 
in the Supplementary Information. The resulting scale for the displacement 
spectral density Sxx is shown on the right side of Fig. [^. Note that the chosen 
convention of x allows directly comparing the readout of the two mechanical 
resonances on this scale. 

To determine the sensitivity with which the displacement spectrum of the 
beam can be read out, we consider the detection noise floor for the measurement 
shown in Fig. which is composed of electronic noise of the photodetector and 
the optical shot noise of the detected light. Their measured combined impreci¬ 
sion (blue datapoints in Fig.|^) is over 7 orders of magnitude smaller than the 
measured signal. 

A general assessment of the sensitivity capabilities of the measurement is 
made by comparing the detection noise imprecision to the (shot noise) impreci¬ 
sion 5'^™P(f]rn) of a resonator read out at the standard quantum limit (SQLjUEH 
The imprecision at the SQL is equal to half of the spectral density of the zero- 
point fluctuations /2. We determine this value from the measured ther¬ 

mal noise spectrum of the lowest-frequency mode via the average phonon occu¬ 
pancy of the mechanical mode ksT and indicate it in Fig.[^ with the red 
dotted line. The optical shot noise of the light impinging on the detector, and 
even the total measurement noise floor, are lower than the imprecision noise at 
the SQL. 

Readout of a nanomechanical resona tor with an imprecision below that at 
the SQL was first achieved in 200^53111 making use of high-quality optical and 
mechanical modes. These high quality factors were instrumental because the 
ability to perform a measurement with SQL-level sensitivity scales, per intra¬ 
cavity photon, with the single-photon cooperativity Co = nV. This shows it 
depends on the photon-phonon coupling strength as well as the optical linewidth 
K. = cJc/Qopt and the mechanical linewidth F = The fact that here we 

achieve a detection noise imprecision below that at the SQL with optical and me¬ 
chanical quality factors of both less than 500 attests to the large optomechanical 
coupling strength, and could have important application in broadband, sensitive 
nanoscale sensors. 

Determining the photon-phonon coupling rate. To quantify the op¬ 
tomechanical interaction strength in the fabricated devices, we model the trans¬ 
duction of thermal displacement fluctuations using equation Q and use it to fit 
a low-power measurement on a structure for various laser detunings. We do this 
by calculating the variance of the optical power fluctuations 5P at the detector 
resulting from displacement fluctuations 5x of a mechanical mode with known 
(thermal) variance. Integrating equation (§ over a single mechanical mode and 
using our expression for the reflection spectrum R(A) (equation 0), yields 

C,p2 2 ksT (A^ex “ cAkcos(</?) - c(A2 - /A) sin((^))^ 

(A2 + kV4)4 

which is independent of the choice of the displacement coordinate x. 


(p2 
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Figure 2 | Pree-space characterisation, (a) Schematic diagram of the free-space 
readout method (PBS: polarizing beamsplitter; H,V: horizontally and vertically polar¬ 
ized light. See Methods for details), (b) Reflection spectrum (red datapoints) and fit 
with a Fano lineshape (black line), (c) Modulation spectrum of the reflected light ob¬ 
tained with the laser frequency on-resonance with the cavity (orange datapoints), and 
a fit of the two mechanical resonances (brown line). The noise floor (blue datapoints) 
was obtained by reflecting the laser light from the unpatterned substrate and matching 
the intensity on the detector. The black dashed line is the sum of the measured elec¬ 
tronic noise (grey datapoints) and the optical shot noise calculated from the intensity 
on the detector (light blue dash-dotted line). The red dotted line shows the peak value 
of Sp ^/2 for the lowest-frequency resonance, which we obtained from the fit of the 
measured thermal spectrum via the relation *S'p^p(Qm)/2 = *S'^p(Qm)/if^m/4/csT. 


The measured variance of the optically modulated signal due to the lowest- 
frequency mechanical mode is shown in Fig. [^. The variance is minimal when 
the derivative of the reflection signal (Fig.|^) vanishes. Interestingly, due to the 
dispersive lineshape the transduction is largest for the laser tuned to resonance. 
The line shown in Fig. is a fit of equation to the data, using only go as 
a free fitting parameter (all other parameters having been determined in inde¬ 
pendent measurements). The corresponding value for gol27r is 11.5 MHz, which 
is an order of magnitude larger than previously reported valueSlElElElHlH 

To compare this photon-phonon coupling rate to the prediction from our 
simulation we estimate the zero-point fluctuations of the structure. Using the 
measured mechanical resonance frequency and the simulated effective mass, we 
obtain x^pf = ^ 0.08 pm for a purely anti-symmetric mode. With 

the simulated frequency response G, this yields a prediction of goj^n 26 MHz. 
To take into account the observed asymmetry of the mechanical mode, we should 
apply a correction factor of 0.76, based on our knowledge of 0 (see Supplemen¬ 
tary Information). This results in an expected value of got^n ~ 20 MHz. We 
attribute the remaining discrepancy to fabrication imperfections, that could re¬ 
sult in a different overlap of the optical and mechanical modes than simulated. 
So in fact, these simulations show that a further increase of go even beyond the 
measured value is possible. 

Optical spring tuning. While we tune the laser frequency across the op¬ 
tical resonance a pronounced shift of the mechanical resonance frequency is 
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Figure 3 | Sensitivity of the displacement measurement, (a) Reflection spec¬ 
trum (red datapoints) and fit with Fano lineshape (black line) for this particular 
nanobeam, (b) Detected optical variance in the fundamental mechanical resonance 
(blue datapoints) with power incident on the sample Pin = 8.5 pW. The datapoints 
were obtained by fitting the fundamental mechanical resonance peak in the measured 
modulation spectra (inset). The signal originates from thermal motion so it varies only 
with the sensitivity of the measurement. The black line shows our model, which uses 
the parameters obtained from the reflection spectrum in a and is fitted to the data to 
determine go — 11.5 MHz. 


observed. In Fig. this is shown for the same structure we studied in Fig. 
This well-known optical spring effect is caused by the radiation pressure force 
being opposed to (aligned with) the mechanical restoring force when the laser 
is detuned below (above) the resonance frequency, changing the effective spring 
constant and therefore the mechanical resonance frequenc}^^. The equation that 
describes this behaviour in the limit of a large cavity linewidth {n ^ Vt^) is 

_ 2 ^in-Pin 

”“^ 0^2 + ^ 2/4 

From equation Q we recognize that the optical spring tuning shown in Fig.[^ 
provides a second, independent way to characterize the photon-phonon coupling 
rate. Figure]^ shows the center frequency of the mechanical resonance extracted 
from the same measurement as the variances in Fig. [^, as well as a fit using 
equation To estimate qq from this fit we need to know /^in, which we cannot 
easily determine as it generally depends on the overlap between the focused 
Gaussian beam and the cavity mode profile. However, we can find bounds for 
/^in by considering the total decay rate n and /^ex = \//^in/^out 5 which were 
determined from the fit to the reflection spectrum. On the one hand we know 
I-C. the collectiou efficiency is at least as efficient as the overlap with 
a Gaussian beam, and on the other hand K[n > i.e. at most half of the 

light escaping from the cavity can be collected because of the vertical symmetry 
of the structure. Gombining these bounds with the fit of the optical spring effect 
yields a range for go between 10 and 13 MHz, in good agreement with the value 
obtained from the analysis of measurement transduction. The fact that the 
spring shift can be fully explained by the radiation pressure force as predicted 
by equation Q shows that forces due to photothermoelastic effectare likely 
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Figure 4 | Optical tuning of the mechanical resonance frequency, (a) Spectro¬ 
gram showing optical tuning of mechanical resonance frequency with power incident 
on the sample Pin = 140 pW. (b) Fitted frequency of the fundamental mechanical 
resonance (blue datapoints) with Pin = 8.5 pW. The black line shows a fit using the 
model for optical spring tuning (equation 0 )- 


insignificant compared to radiation pressure. 

Nonlinear transduction. As a consequence of the large photon-phonon 
coupling rate, the thermal motion of the nanobeam induces frequency changes 
that are appreciable with respect to the linewidth of the cavity, which results in 
nonlinear transduction. This generates spurious signals at integer multiples of, 
and combinations of, the strongest modulation frequencies. Detection of such 
signals at multiples of the mechani cal reson ance frequency resulting from ther¬ 
mal motion was reported previousl^EHMIMl and compared to quadratic optome¬ 
chanical couplin^^. 

Figureshows a transduced spectrum where we identify 15 peaks as integer 
multiples and combinations of the two fundamental mechanical resonances at 
1.4 MHz (“A”) and 2.0 MHz (“B”): = IjA ± /cB|, with j, /c G {0,1, 2,... }. 

Peaks corresponding to different order {j + k) have a different detuning depen¬ 
dence, but all peaks with the same order differ only by a constant factor. To 
illustrate the detuning dependence of the higher-order peaks, we plot the vari¬ 
ance of the peaks jA for j = {1, 2, 3,4} in Fig. [^. The detected height of the 
higher-order peaks can be predicted by a Taylor expansion of the amount of light 
in the cavity around the average detunin^^ (see Supplementary Information), 
the result of which is shown in Fig. [^. 

Note that the higher-order peaks in this calculation were not fitted to the 
data, but follow from the value of go we obtained by fitting the first-order 
peak, as shown in Fig. The measured nonlinear sidebands are larger than 
expected (corresponding to a suggested increase of go of about 60%). The origin 
of this discrepancy is unknown. Possible explanations include higher-order op¬ 
tomechanical couplin^^or mechanical nonlinear it ie^^. However, the symmetry 
and shape of the curves match the experimental data, which confirms that the 
detuning dependence corresponds to the successive derivatives of the reflection 
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Figure 5 | Nonlinear transduction, (a) Measured mechanical spectrum of the 
sliced nanobeam (purple datapoints), using near-resonant light with power incident 
on the sample Pin = 8.5 pW. Three peaks A, B and C are mechanical resonances of 
the nanobeam, all other visible peaks correspond to integer multiples and combinations 
of frequencies A and B. The electronic noise floor is shown with grey datapoints. (b) 
Areas under the peaks corresponding to integer multiples of frequency A, obtained by 
htting the peaks in the spectrum, (c) Calculated variance of the reflected signal, using 
the experimentally obtained parameters. 


spectrum (Fig. [^). 

Discussion. The free-space readout method we employ provides an easy 
and robust way of coupling light to the cavity. We have intentionally engineered 
the cavity defect such that it has a significant dipole momenfP^, allowing cou¬ 
pling to free space at an appreciable rate. This makes it unnecessary to create 
an explicit loss channel for coupling, e.g. in the form of a grating or feeding 
waveguide. Moreover, the Fano-shape of the reflection spectrum allows direct 
transduction of motion to optical amplitude modulation for a laser tuned to 
the cavity resonance (where dynamical radiation pressure backaction is zero), 
without more complicated interferometric schemes. As a result of the efficient 
coupling to free space, the bandwidth of the cavity is large (0.5 THz), which 
is appealing in the context of applications that require frequency matching of 
multiple systems: together with the small system footprint, it could assist the 
integration of such optomechanical transducers in sensor arrayP^ or effective 
optomechanical metamaterialJ^. 

Of course, for applications that benefit from enhanced measurement sensi¬ 
tivity such as measurement-based control of the mechanical quantum state, it 
could be worthwhile to reali ze a higher optical quality factor by introducing 
tapering along the nanobeanlllESl To simultaneously allow efficient free-space 
coupling would in such a case require special attention, in the form of tailoring 
the spatial mode profile of the cavity radiation. This could be especially im¬ 
portant for effects that depend on the intracavity photon number, such as the 
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demonstrated optical spring effect, as the rates and /^out will differ. Further 
quantification of their individual magnitudes (e.g. through systematic variation 
of incident and detected mode profiles) will thus be valuable. 

Likewise, we expect that the mechanical quality factor for the nanobeams we 
employ can be improved with suitable design principles and optimization of the 
fabrication process. Indeed, measurements on similar-sized silicon nanobeams 
and cantilevers suggest that quality factors in the range of 10^ to 10^ should 
be possible at room and cryogenic temperature, respect ivel^HEH Nonetheless, 
we point out that because of our large coupling rate, even with the current 
modest values of both optical and mechanical quality factors the single-photon 
cooperativity in this structure reaches Cq = 0.16. The value of this quantity, 
which compares optomechanical coupling strength and dissipation, and is for 
example a measure for the capability of the system to perform measurements 
at the SQL, is on par with many recently reported systems with much higher 
quality factor^. 

In conclusion, we demonstrated an optomechanical device with a large photon- 
phonon coupling rate ^o/27r = 11.5 MHz, and used it to demonstrate sensitive 
measurement of nanomechanical motion and pronounced optical tuning of the 
mechanical resonance frequency. It is interesting to note that the regime of 
large coupling rate and modest optical linewidth is beneficial in the context of 
achieving strong mechanical tuning, as parametric instability is suppressed. We 
revealed that the working mechanism relies on an optical mode with a subwave¬ 
length mode area. We predict this approach can be extended to yield even larger 
coupling rates, or to be applied to modes with higher mechanical frequencies. 
In the current device the photon-phonon coupling rate go exceeds the mechan¬ 
ical res onance frequency (fm, which is one of the requirements for ultrastrong 
couplin^lElHH With further improvements in both the coupling rate and the 
optical quality factor, the present approach might provide a route to simultane¬ 
ously reach go > and go ^ It will be interesting to explore to what extent 
this regime can be used to exploit nonlinear optomechanical interactions at the 
single-photon level. 
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Methods 

Numerical simulation. All numerical eigenmode simulations were performed 
using finite-element software COMSOL Multiphysics. In mechanical simulations, 
the connection between the substrate and the support pads was modeled as a 
fixed boundary, while all other boundaries were kept free. To find the guided 
modes of the photonic crystal nanobeam, a unit cell was simulated with Floquet 
boundary conditions along the propagation direction and in the other directions 
perfect electric conductors at several micrometers distance from the structure. 
Finally, to simulate the cavity mode, a full nanobeam including support pad was 
modeled with perfectly matched layers on all sides, again at several micrometers 
distance from the beam. 
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Fabrication. The structures were fabricated from a silicon-on-insulator 
wafer with a device layer thickness of 200 nm, and a buried oxide layer of 1 pm 
thick. A resist layer of hydrogen silsesquioxane (HSQ) with a thickness of 80 nm 
was spincoated on top and patterned using electrons accelerated with 30 kV. 
The resist was developed using TMAH and then the pattern was transferred 
to the silicon layer using a reactive-ion etch process with SF 6 /O 2 gases, opti¬ 
mized for anisotropy and selectivity. To release the structures, the oxide layer 
was dissolved in a 20% HF solution. After this step the structure was dried 
with a critical point dryer to prevent the sliced beams being pulled together 
during the drying process. The suspension of the nanobeams from their support 
pads was designed to allow some relief of compressive stress along the beam. 
The compressive stress is present in most free-standing structures created from 
SOPl, but it has a large effect for our structures because of their low stiffness. 

Pree-space setup. The laser beam (New Focus Velocity 6725) was focused 
on the sample by an aspheric lens with a numerical aperture of 0.6. A polarizing 
beamsplitter provided a cross-polarized detection scheme, where any light that 
was directly reflected was rejected and only light that coupled to the sample, 
placed at 45°, was transmitted to the detector. Both the lens and the sample 
were in a vacuum chamber to reduce mechanical damping by air molecules. All 
experimental results shown were performed with a pressure of about 4 mbar, 
except the spectra in Fig.[^ where the pressure was lower than 10“^ mbar. At a 
pressure of 4 mbar, the mechanical quality factor was lowered to approximately 
200 . 

Analysis of modulated reflection signals. We detected the reflection 
signal using a low-noise InGaAs-based photoreceiver (Femto HCA-S) and an¬ 
alyzed it using an electronic spectrum analyzer (Agilent MXA). We fitted the 
peaks in the modulation spectra using a Lorentzian convolved with a Gaussian 
distribution, also called a Voigt lineshape. The Gaussian contribution accounted 
for the resolution bandwidth of the spectrum analyser, as well as for frequency- 
noise broadening at relatively large optical input power Pin- With high Pin, small 
fluctuations in incoupling efficiency or laser intensity thermally shifted the cav¬ 
ity resonance, which resulted in frequency noise via optical spring tuning of the 
mechanical resonance. 
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Supplementary Information 
Model for transduction 

We derive the effect of a small frequency modulation of the intracavity field at 
the output of a cavity, which allows us to set up a model that describes the trans¬ 
duction of thermal motion to optical intensity modulation. This derivation fol¬ 
lows the calculation shown by Gorodetsky, Schliesser, Anetsberger, Deleglise & 
Kippenberg , with the important difference that we include the non-resonant 
contribution in the reflection spectrum that leads to the Fano lineshape we 
observe. 

We start from the equations describing the behaviour of the optomechanical 
system 


a = (i (A - Gx{t)) - hi/2)a{t) + ^/K^Sin{t), 
x{t) + Tmx{t) + D^x(t) = -hG\af, 


(SI) 


where a is the internal field in the cavity, S[n is the input field related to the 
input power Pin = huj\s[n\‘^, A is the detuning of the input light from the cavity 
resonance cdc, ^ is the cavity decay rate, /^in is the coupling rate to the input 
channel and G = duoddx is the optomechanical frequency response. The fre¬ 
quency Dm and damping rate F^ of the mechanical resonator are influenced by 
the number of photons in the cavity |dp, an effect we neglect in the following 
by assuming a low input power. This simplification is motivated by the fact 
that we seek to predict the amplitude of the mechanically-induced light modu¬ 
lation, not the frequency of such modulations. Moreover, dynamical backaction 
affecting mechanical linewidth is small in devices that have large /^/Dm- 

We consider a small harmonic oscillation of the mechanical resonator x{t) = 
xq cos(Dmt), which causes a modulation of the cavity frequency with amplitude 
xqG, or a modulation of the optical intracavity phase with amplitude xoG/Dm- 
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If the modulation is small {xqG /^), this yields 
— <5in a/^ in ^(0) 

X (^1 - ^-xoG , ^g2) 

^ -i{A + n) + K/2' 

The cavity is coupled to the output 5out with the coupling rate /^out- The equa¬ 
tion for the output field reads 

•^out — CC^SiYi \/ ^out^a: 1 

where the first term is caused by the nonresonant scattering from the input to 
the output with amplitude c and phase (f. 

In the experiment, we measure the intensity l^outP and feed it to a spectrum 
analyser, which yields the single-sided spectrum of the signal. Since we are 
interested in the strength of the spectral component at the mechanical oscillation 
frequency we highlight the time dependence here: 

|SoutP(i) = c^kinp + Kout\ax{t)\'^ “ (e*‘^Sin< (*) + s*^a^{t)). (S4) 


The first term is constant so it does not contribute to a signal at Substi¬ 
tuting ax into the other two terms, and discarding any terms not oscillating at 
yields 


|SoutP(i) 






X m? - £(-f^m))] + C.C.) 


-c([ie*“"'‘ (e*^£*(0)£*(^^m) -e-*‘^£(0)£(-^^m))] + c.c.) 

(S5) 

This expression contains the modulation amplitude. In our experiment, we com¬ 
pare the variance of the modulation to the known variance of the mechanical 
thermal motion (x^)th = /TiVt^. Therefore we calculate the variance 

of Pout = ^<^o|5outP due to the modulation at and —which will both 

contribute to the signal at in the single-sided spectrum. We can write 

Pout I = which leads to (|PoutP)om = 2|x4p. After some 




algebra, we arrive at 


{\Pont?)a^ = 


l^in ^out) 


(A^ -h /a:2/4)2((A + + /a:2/4)((A - (2^)^ + n^/A) 

A i^in^out 2Acy^COS ip -|- (A ^ / d) siu ip^ 

coA^{ip) 

+ (A^ — A^/a:^/2 + /lA) sin^((/?) 

— 2 cos((/?) sm{(p){—A^K + Ak^/4 + AkCI‘^/2) 


(S6) 
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If we evaluate this expression in the bad-cavity limit ^ we find it 
is directly related to the derivative of the Fano lineshape dR/dA: 


■ out I 


— 

/Om — 


2 xlG'^P^^KinKout (Ay'KinKout “ cAkCOs{(p) - c{A'^ - /P) sin((yj))" 

(A^ + /^^/4)^ 


2 0*^ ^in I 


2 


(S7) 

We note that imperfect transmission of the optics between the sample and 
the detector scales the detected signal in the same way as the input power Pin, 
and will enter the equations in the same way. Finally, we substitute the variance 
due to the modulation amplitude xq by the variance of the mechanical motion: 
XqG‘^ 2(x^)thG^ = which leads to equationj^in the main text. 


Nonlinear transduction 


The previous section started from the assumption that the frequency modula¬ 
tion dujc = GSx is small with respect to the cavity linewidth, SuJc ^ /^, and 
considered only the resulting linear transduction at the modulation frequency 
In this section we show that the first signature of large 6uJc is the ap¬ 
pearance of nonlinear transduction, which produces a signal at multiples of the 
modulation frequency For the second-order transduction, this was shown 
by Doolin et al. [^, where also a quadratic optomechanical coupling was taken 
into account. Here we derive the result for any higher-order terms of nonlinear 
transduction. 

In the non-resolved sideband regime {n ^ the optical fields in the cavity 
reach a steady state much faster than the timescale of mechanical motion. The 
intracavity amplitude can then be written as 


a{t) = 


/^in^in 


—i(A — 5ujc{t)) + n/2'‘ 
which combined with equation (|S3|) yields 


|<5out I 


2 I d/^in/^out 

= c + -^- 


e*^(l - m)-h e 


1 


(S8) 


(S9) 


Here we defined u = 2{A — SuJc{t))/K, which implies u is detuning- and time- 
dependent. We now summarize equation (S9) as R'{u) and find the Taylor ex¬ 
pansion for small 6uJc around uq = 2A/n: 


R\u) = R'{uo)- 


2SuJc dR'{uo) 
n du 


{-2SuJc/k)^ d^R'{uo) 


n\ 


du^ 


(SIO) 


where the last term depicts the nth order in the Taylor expansion. We take 
a harmonic modulation of the cavity frequency 6uJc{t) = A cos fit. To leading 
order, Suj^ ^ A'^ cos{nfit)/2^~^. This means that each successive term in the 
Taylor expansion in equation (SIO) gives the amplitude of a term at different 
frequency. 

In the optomechanical system, the variance of the frequency modulations at 
the mechanical frequency fl^ is given by {Suj‘^) = G^(x^)th = 2gQkBT/Tifl^. 
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Therefore we get the same variance in R if we set A = 2go^/kBT/M^m‘ We note 
that ^ ^ which for the variance of the signal at leads 


to 


(P2) 




nQri 



2{glkBT/nn^r (d^R{A)Y 

nP V dA^ J 


(Sll) 


For n = 1, the result of equation in the main text is again reproduced. The 
result of a calculation of the variances of n = 1-4, for the parameters used in 
the experiment, is shown in the main text in Fig. [^. 


Waveguide modes in the sliced nanobeam 

In this section we discuss the waveguide modes in the periodic region of the 
sliced nanobeam in more detail. The free-standing silicon nanobeam acts as a 
waveguide, which can guide light via total internal reflection. In this waveguide, 
the elliptical holes form a photonic crystal that opens a bandgap for modes 
with transverse electric (TE)-hke symmetry (Supplementary Figure [S^). This 
is not a full bandgap, since TM-like waveguide modes exist in the gap region. 
If the symmetry of the structure is broken by fabrication imperfections, light 
in the bandgap region for TE-like modes can scatter to the TM-like modes and 
propagate along the nanobeam. Eor this reason, this is sometimes referred to as 
a quasi-bandgap. 

The guided mode at the lower edge of the band gap has the largest concen¬ 
tration of energy in the nanoscale gap in the middle of the beam. The fact that 
a significant portion of this mode’s energy is located in vacuum increases its fre¬ 
quency in comparison to a non-sliced nanobeam, which reduces the frequency 
width of the bandgap. To ensure maximum mirror strength, the transverse size 
of the holes is made as large as possible. The elliptical hole shape is as such 
important to realize a strong bandgap, in addition to providing favorable me¬ 
chanical properties as mentioned in the main text. 

To create optical cavity modes that are derived from the lower band edge, 
the defect is a local decrease in distance between two elliptical holes, which 
decreases the local effective refractive index and creates defect states in the 
bandgap region. Supplementary Eigure shows the first two cavity modes 
created in this way in the sliced nanobeam. Note that the higher-order cavity 
modes have a lower frequency, since they are less confined near the defect, so that 
the frequency is closer to that of the waveguide mode in the periodic structure. 


Overlap of optical and mechanical mode profiles 

The frequency shift of the optical resonance of the sliced nanobeam due to me¬ 
chanical motion depends on the overlap between the optical and mechanical 
mode profiles. We simulated the frequency shift of the optical resonance of the 
sliced nanobeam due to a uniform mechanical shift of I nm, as shown in the 
main text in Eig. To estimate the influence of the finite extent of the mode 
profiles on the response, we extract these profiles from the numerical simulation. 
Supplementary Eigure [S^ shows the normalized displacement profile of the me¬ 
chanical resonance and the normalized electromagnetic energy density profile of 
the optical cavity mode, as a function of the position along the nanobeam. 


18 












k\\[nl a) 


Figure SI | Optical modes in the sliced nanobeam, (a) Simulated disper¬ 
sion diagram showing the TE-like waveguide modes in the periodic part of the sliced 
nanobeam. A bandgap is opened by the periodic structure of elliptical holes around a 
frequency of 200 THz. (b) Simulated transverse electric field prohles of the first two 
cavity modes of the structure. The defect that is responsible for the creation of these 
modes is a slightly smaller distance between the two holes in the center of the beam. 


The mechanical mode profile closely resembles the mode profile of the fun¬ 
damental mode of a doubly-clamped beam. From the simulated displacement 
profile we calculate the effective mass for this purely antisymmetric motion to 
be rueff ~ 0.39m, where m is the total mass of the sliced nanobeam, which is 
indeed very close to the value obtained from the analytical displacement profile 
of a uniform doubly-clamped beam[^. 

The optical mode profile clearly shows the localization of the energy density 
in the small gaps between the silicon ‘teeth’ of the structure. The simple defect 
we introduce in the center of the beam localizes the optical cavity mode there, 
while the field decays exponentially away from the defect, where the optical 
frequency lies inside the photonic quasi-bandgap. 

We compute the correction on the frequency shift due to the finite extent of 
the modes from the overlap integral between these two mode profiles, and find 
a factor of 0.90 with respect to the frequency shift for a uniform displacement 
of the beam. 

The higher-order cavity modes created by the defect are less strongly con¬ 
fined along the length of the beam, as shown in Supplementary Figure [STj p. 
Since the fundamental mechanical resonance has the largest displacement in 
the center of the beam, the higher-order cavity modes are less sensitive to this 
motion than the fundamental optical resonance. 


Mechanical mode coupling 


The mechanical modes observed in a two-beam system are combinations of the 
motion of the individual beams, with in the perfectly symmetrical case fully 
in-phase or out-of-phase motion. Here we derive the consequences of imperfect 
symmetry for the ratio of scattered power between the two modes [^. 

Harmonic motion of the two beams at a certain frequency 12 can be described 
as: 


Xi{t) 


'tpl 

_ X2{t) 




cos fit = cos 12t, where -0 = 


^2 


(S12) 


Thus, jjji and 7/^2 are the amplitudes of the oscillatory motion of the two indi- 
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Figure S2 | Normalized resonant mode profiles of the fundamental optical and me¬ 
chanical resonance of the sliced nanobeam. The red dotted line shows the displacement 
along the length of the beam, while the blue solid line represents the local energy den¬ 
sity. 


vidual beams, such that their variance is 2 ) = 2 - 

The optical response of the system is determined by the change in the dis¬ 
tance between the beams: d = xi— X 2 . Here we define x as x = d/2, which leads 
to the same value of G = duoddx for both mechanical modes. Note that this 
choice of x corresponds to the lab-frame displacement of the two beams if they 
move in antiphase. The variance of x due to harmonic motion described by ^ is 
then: 

{x^)^ = i (V’l + ■*/’2 - 2 '^iV’2) • (S13) 

The state vectors of the two normal modes can be written without loss of 
generality as 


= Aa 


cos6> 


sinO 

sin^ 

II 

— cosO 


(S14) 


For both of these modes, we can calculate the variance of x, denoted as {x‘^)a 
and {x^)i 3 , respectively: 


(cos^ 0 + sin^ d + 2 sin 0 cos 0) 

8 

= ^(1 +sin 261), (S15) 

8 

{x^)l3 = -sin26») 


For the beams undergoing thermal motion, the variance is given by the 
equipartition theorem: 


2 \ ^bT , 2 \ ^bT 
a = - {X )/3 = - 


(SI6) 


where rrio, and are the effective mass of these modes. As shown in the 
previous section, for the differential mode the simulated effective mass with 
respect to the displacement coordinate x is rueff = 0.39m, with m the total mass 


of the sliced nanobeam. Evaluating equations (S15) and (S16) for a differential 
mode {0 = 7r/4) yields = AkBTj and similarly for A|. Substituting 
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this back into equation ( |S15[ ), we arrive at 



(S17) 


We note that thermal variance is related to the zero-point fluctuations as 



(S18) 


where f3 and + respectively — is chosen as the sign for the term sin26>. 

Finally, we measure (P^) and wish to relate this to a displacement vari¬ 
ance {x^). We calculated rueff using the simulated mode profile and assume the 
thermal bath temperature of the mechanical modes T is equal to the lab tem¬ 
perature, which leaves only 0 and the transduction factor dPjdx unknown. By 
measuring the area of both peaks (P^)a and (P^)/ 3 , we resolve the remaining 
ambiguity, allowing us to calibrate the displacement spectrum. 

Influence of compressive stress and experimental 
disorder 

The simulation of the ideal structure shown in the main text predicts the res¬ 
onance frequency of the fundamental in-plane resonance to be 6 MHz, and ad¬ 
ditionally the frequency difference between the anti-symmetric and symmetric 
mode is negligible. We expect the frequency of the out-of-plane resonances to 
be larger than that of the fundamental in-plane mode, as the narrowest part of 
the half-beams (80 nm) is smaller than the Si slab thickness (200 nm). 

In our experiment we find significantly smaller values of 2.6 and 3.2 MHz 
for one structure, and 1.4 and 2.0 MHz for a second structure. The differences 
between the frequencies of the two in-plane resonances are also larger than ex¬ 
pected. To find the origin of this effect we measure the experimentally observed 
disorder produced by the fabrication process of the second structure using a 
scanning electron microscope. We then introduce these dimensions into the sim¬ 
ulation, which results in slightly different eigenfrequencies. However, the result¬ 
ing values are still near 6 MHz and additionally the two fundamental modes 
still have only a small frequency difference. 

Finally, we add an extra step to the simulation to include a static com¬ 
pressive stress, which can result from the stress in the buried oxide layer of 
the silicon-on-insulator layer structure we use[^. We define an initial displace¬ 
ment for one of the support pads along the length of the beam, which creates 
the stress. Increasing the amount of displacement, and therefore the stress, de¬ 
creases the eigenfrequencies predicted by the simulation. This is in accordance 
with the theoretical expectation]^. At the same time the difference between the 
eigenfrequencies of the two fundamental modes increases. 

Supplementary Figure shows the displacement profiles that result from 
such a simulation. At a displacement value of 10 nm (0.09% of the length of 
the beam), the resonances occur at frequencies very close to the experimentally 
measured values, at 1.5 and 1.9 MHz. 
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Figure S3 | Simulation including disorder and compressive stress. Simulated 
mechanical displacement profiles of the two fundamental in-plane resonances. The 
dimensions of the beam were matched to the realized dimensions using measurements 
with a scanning electron microscope, including differences in hole and gap size along the 
beam. In addition, a compressive stress was introduced in the simulation by displacing 
one of the support pads by 10 nm along the direction of the beam. 


Influence of optical input power 

In the measurements shown in the main text, we use up to 370 pW of optical 
power incident on the nanobeam. Using the parameters of our fit to the re¬ 
flection spectrum, we estimate that this results in an intracavity intensity that 
corresponds to a maximum of ~ 1000 photons. For this intensity in the cavity, 
we do not expect an increase in the cavity temperature of more than a few 
Kelvin. As a first confirmation of this, we see thermal shifts of the cavity reso¬ 
nance frequency of less than 1 nm, which corresponds to a temperature increase 
of less than 10%. 

To check the assumption more thoroughly we measured the signal strength 
of the two fundamental mechanical frequencies as a function of power incident 
on the structure. Supplementary Figure [S4| shows the result for resonant light, 
both from the zero-detuning point in swept measurements and from individual 
measurements where the detuning was set to zero by minimizing the optically 
induced shift of the mechanical frequency. Comparison with the line shows that 
the datapoints closely follow the expected quadratic dependence on input power 
for both peaks in the spectrum. The largest source of uncertainty in this mea¬ 
surement are differences in coupling efficiency between the incoming laser beam 
and the cavity due to small changes in alignments, which can influence the signal 
strength between measurements. 
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